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Trapped-ion Lissajous trajectories 
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Here we present a protocol for generating Lissajous curves with a trapped ion by engineering 
Rashba- and the Dresselhaus-type spin-orbit interactions in a Paul trap. The unique anisotropic 
Rashba a^,, Uy and Dresselhaus I3x, Py couplings afforded by our setup also enables us to obtain 
an “unusual” Zitterbewegung, i.e., the semiconductor analog of the relativistic trembling motion 
of electrons, with cycloidal trajectories in the absence of magnetic fields. We have also introduced 
bounded SO interactions, confined to an upper-bound vibrational subspace of the Fock states, as an 
additional mechanism to manipulate the Lissajous motion of the trapped ion. Finally, we accounted 
for dissipative effects on the vibrational degrees of freedom of the ion and find that the Lissajous 
trajectories are still robust and well defined for realistic parameters. 
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Increasing interest in quantum simulations and con¬ 
trollable trapped ion systems —used through the last two 
decades as a staging platform to investigate fundamen¬ 
tal quantum phenomena tl| and to implement quantum 
information processing [2| — have led to the emulation 
of the quantum relativistic wave equation. Apart from 
trapped ions Si, a number of distinct setups have been 
used to simulate zitterbewegung in a variety of physical 
systems such as quantum wells Q , photonic crystals i , 
Bose-Einstein condensates i and ultracold atoms Q ■ It 
has been demonstrated that Bose-Einstein condensates 
can provide an analog to sonic black holes Si ^.nd the 
Higgs boson 0. _ 

In trapped-ion system Gerritsma et al. [HI reported 
a proof-of-principle quantum simulation of the one¬ 
dimensional Dirac equation within a trapped-ion exper¬ 
iment. Measuring the time-evolving position of a single 
trapped ion set to mathematically simulate a free rela¬ 
tivistic quantum particle, they demonstrated the Dirac 
Zitterbewegung, as anticipated in Ref. S, distinct 
initial superpositions of positive- and negative-energy 
spinor-like states. The momentum-spin coupling, which 
naturally takes place in the relativistic quantum regime, 
was simulated by the laser-induced coupling between the 
vibrational and the electronic states of the ion. 

Motivated by the exciting results in Refs. iim , here 
we consider a single two-level atom in a two-dimensional 
Paul trap with four electrodes, Eig.l, and show how to 
obtain spin-orbit interactions of the Rashba and Dres¬ 
selhaus types, ubiquitous in quantum spintronics and in 
emerging areas of condensed matter physics such as topo¬ 
logical insulators and Majorana fermions. Our setup en¬ 
ables us to engineer Hamiltonians with anisotropic SO 
couplings such as 

TT PpE 0/jj Px , Py /I \ 

Hso = ^Cry,Py - —ayPx + ^CTxPx + ^CFyPy, (1) 

which generalize the canonical form of the SO interac¬ 
tions in semiconductor quantum wells to include distinct 



FIG. 1: The schematic of the experimental Paul trap setup 
illustrating the ion performing a Lissajous curve. 


Rashba ax, ay and Dresselhaus j3x, Py couplings. In (HD 
p and (T denote the momentum and pseudo-spin opera¬ 
tors, respectively, of the trapped ion. Interestingly, the 
anisotropic SO couplings can give rise to physical phe¬ 
nomena not possible in the condensed matter environ¬ 
ment, e.g., ionic motion with trajectories following Lis¬ 
sajous curves, Eig.l. In Fig. 2 we show the whole set 
of Lissajous curves followed by the trapped ion, which 
we discuss in detail below. In addition, we can also ob¬ 
tain the “unusual” Zitterbewegung, a trembling relativis¬ 
tic motion, wit h cy cloidal trajectories in the absence of 
magnetic fields [1^. We also demonstrate how to engi¬ 
neer bounded-SO Hamiltonians, confined to an upper- 
bound vibrational subspace of the Fock states —leading 
to distincts periodic motions— apart from analysing the 
effects of dissipation over trajectories. 

Engineering the anisotropic Rashba interaction. We 
consider a two-level ion acted upon by two laser beams 
perpendicular to one another in x and y directions, each 
providing red and blue sidebands simultaneously via an 
electro-optical modulator [l^. From the four generated 
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frequencies and relative phases ipi, two are tuned to 
the first-blue sidebands, each coming from one of the orig¬ 
inal laser beams. The other two are tuned to the first-red 
sidebands, leading to the Hamiltonian {h= 1) 

+^yr]y -I- cr+ -|- H.c.] , (2) 

where rja is the Lamb-Dicke parameter and the 
Rabi frequency for the coupling between the electronic 
(ground \g) and excited |e)) states of the ion and its two- 
dimensional motional degrees of freedom, described by 
the annihilation as and creation a| operators with respect 
to each direction s = x,y. Moreover, cr+ = |e) {g\ and 
cr_ = cr(^ are the raising and lowering operators for the 
two-level system. By adjusting (pi — 0,<p2 = tt, (pa = 7r/2, 
and p 4 = 3ttI2, we obtain, in the interaction picture, the 
anisotropic 2D Rashba SO interaction 


which, together with the Rashba coupling, is a crucial in¬ 
gredient to manipulate the electron spin in quantum spin- 
tronics and, more recently, to obtain exotic topological 
phases of matter. We can also consider four laser beams 
simultaneously and hence generate a Rashba-Dresselhaus 
Hamiltonian as shown in Eq. o, with as = AsTy^Ds and 
Ps = As77sDs, enabling us to simulate the canonical SO 
interaction of quantum spintronic in our setup. 

At this point it is important to stress that a four-level 
ion, as in Ref. Q, can also be used to engineer the in¬ 
teractions of the forms [Eqs. ([T])- ([B])]. In this case we 
have to consider two additional states \g') and |e') having 
the same angular momentum as jg) and |e), respectively, 
but distinct energies and magnetic momenta. For in¬ 
stance, by exciting only the crossed transitions g o e' 
and e -O' g', we can exactly simulate the mathematical 
structure of the Hamiltonian used to obtain the unusual 
zitterbewegung in Ref. 12|. In this case, the isotropic 
Rashba-type interaction ([3]) is modified to 


H - {^/S.xTjx^xVx^y AyT/yHyPyCTa,) , (3) 

where Ps = i (a| — a^) /2As, with As = Ij^2mvs^ ^ 
being the ion mass and Vs the trap frequency in the s 
direction. When we adjust the laser parameters so that 
Ax = Ay, r]x = gy, and fix = Dy, the above derived 
Hamiltonian reduces exactly to the usual Rashba inter¬ 
action i? = 7 {px<Xy — Py<Jx) in quantum wells, where 

7 = Agfl. 

Interestingly, we find that a 3D trap with an additional 
laser beam and an appropriate set of phase constants pi 
enables us to engineer a 3D Rashba-like term 

H = yf. {p X a), (4) 

where r sets the position of the ion in relation to the 
trap center. Unlike the 2D and 3D cases in Eqs. and 
o, the ID Rashba-like interaction H = ^Px<Jy does not 
lead to the unusual Zitterbewegung with cycloidal orbits. 
However, when one of the original laser beams leading 
to Eq. ([2]) is adjusted to produce a carrier interaction 
(instead of the simultaneous red and blue sidebands), we 
end up with the ID interaction 

H = ^Pxay + ^(Jz. ( 5 ) 

Here, as in the Dirac equation, the energy gap provided 
by the Stark shift term Duz (D being also a Rabi fre¬ 
quency), is a key ingredient for producing Zitterbewe¬ 
gung. In addition, an appropriate choice of parameters 
of the laser producing simultaneous blue and red side¬ 
bands leads straightforwardly to a Hamiltonian similar 
to that in Eq. H = ryHAj, Px^x + Dctz, which simu¬ 
lates the ID Dirac equation Q. 

Anisotropic Dresselhaus term. We can also simulate 
the linear Dresselhaus SO interaction by adjusting pi = 
—p 2 = tt/2 and p^ = —pi = 37r/2 in Eq. ([3]); we find 

H = J {pxCTx - PyO-y) , ( 6 ) 


i? = da; ® r {PxO-y - PyCTx) , 


( 7 ) 


where E = AryD and dx is the Pauli matrix in subspace 
iW) (W)}- All other Hamiltonians [Eqs. (H]) - ®] are 
also modified via a tensor products with dx ■ For simplic¬ 
ity, we restrict our analysis to a two-level system which, 
as shown below, gives rise to the same dynamics for the 
unusual Zitterbewegung as given in Ref. [l^ . 

‘'^Bounded” spin orbit interaction. Before analyzing the 
Lissajous curves in Fig. 2, we would like to address a 
case, with no analog in the solid state physics, in which an 
upper-bound interaction 1^ between the internal and the 
vibrational degrees of freedom of the ion is engineered. 
Following the steps outlined in Ref. 1^, we can obtain 


Hk=±i = xiv) [Mv)^± + > (8) 

where xiv) = V ~ ^^/2) and A{g) = 

[l — ?7^ala/2] a. Expanding the operators A and 
in the Fock space basis and adjusting the Lamb-Dicke 
parameter to = 2/N [ry^ =2/ {N — 1) with A G Z] 
such that Al'(ry) \N) = 0 [A(ry) \N) = 0], we find that for 
an initial vribational state prepared within the upper- 
bound subspace ranging from |0) to \N), the Hamiltonian 
Hk=±i becomes 

N-l 

= '^Xn{\n) {n-\-l\a±-\- H.C.). (9) 

n—0 

where Xn = \/n + 1 (l — g^nl2') x{g)- By using the 
Hamiltonian (IT|) instead of ([3]) we can derive all the above 
spin-orbit interactions but restricted to within the upper- 
bound subspace. In particular, for the case of the upper- 
bound Rashba-Dresselhaus interaction, with the same ad¬ 
justment of the laser parameters as above we obtain 

H = axPx^Uy - aypfux + PxPffTx - PyPy^dy, (10) 
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where pf = i [AJ(77) - As{f])] /2As, as = XnV^s, and 
Ps = Xnil^s- should become clearer below, when dis¬ 
cussing the trajectories, the bounded Hamiltonian (fTUl) is 
able to generate Lissajous trajectories other than those 
derived from the usual Rashba-Dresselhaus interaction. 
It also enables Lissajous curves from an initial coherent 
state, which does not occur with interaction ([T|). We fi¬ 
nally note that for each upper-bound parameter N we 
get distincts trajectories. 

Cycloidal Zitterbewegung. Now let us analyze the zit- 
terbewegung with the (anisotropic) Rashba-like coupling 
in Eq. ©, from which we obtain the time-dependent 
components of position operator 

s{t) = s{0)-^^t-iH (^ccrr+ - l) , (H) 


with r s = x,y. We note that Eq. (fTTl) exhibits 
the same time dependence as that coming from the 
Dirac zitterbewegung, with the trembling motion aris¬ 
ing from the third term on the right side of the equal¬ 
ity. Hereafter we consider the initial state |'i/'(0)) = 
f([^exp [—(ps — pos)/2/rJ 0 |p), with both momentum 
eigenstates peaked around pos (ps being the width of the 
distribution) and |p) standing for the ionic internal state, 
we thus compute the mean value (s(t)) = Tr [sp(t)]: 


(s-(r)) = (s-(0)) + (-l)^+''“s^- (a,)r 
epr 


+ 


2^ L 

epr 


T - 1 


i^z) 


sin 


2^3/2 L 
X {{<7x)Px +£(0-y)Py) 


( 12 ) 


This ensures that the cycloidal trajectories without mag¬ 
netic fields, derived in Ref. ([13) by considering the ef¬ 
fects of an interband SO term, can also be obtained here 
from the anisotropic Dresselhaus interaction. However, 
when the Rashba- and the Dresselhaus-type interactions 
come together as in Eq. o, we obtain for the isotropic 
case {ax = ay and Px = jdy) 


(s(r)) = (s(0))+^(-l)^+'-((a.)-«(a,)) 

- Ks>C~'^ (Pr) {aP) [1 - cos {}<t)] 

— Ks>c~^ {pr) * [sin {>ct) — r] 


(14) 


where we have defined the dimensionless cou¬ 
pling K = a/p, and the parameters Ks = 

2 -I-(—1)^“'''^^^^ , = 2k[k^'^+ 2{pxPy)], and 

A = {{Kpy +px) cTy + {npx +py) cTx)- Equation ((T3]) 
shows that with both spin-orbit couplings acting to¬ 
gether with different strengths (k ^ 1), we still obtain 
similar trajectories to those coming from each coupling 
acting separately. Now, when considering matched 
coupling strenghts (k = 1, i.e., a = P), two special 
situations arise when the electronic levels are prepared 
as eigenstates of az or as- In the first case we do not get 
Zitterbewegung as expected, since an eigenstate of Oz 
does not simulate the superposition between positive and 
negative energy states required for the Zitterbewegung. 
However, we do obtain an interesting effect: we lock the 
motion of the particle in the s direction and obtain an 
oscillatory harmonic motion in the r direction: 

(s(r)) = (s(0)) - {pr ), (15a) 

(f(r)) = (r(0)) - Kr>c~‘^ ips) [1 - cos (xt)] . (15b) 


where we have defined the dimensionless position s = 
s/A, momentum ps = 2Aps and time r = rjidt, with 
^ = \/(Ps) + iPl) and e = A^gx^x/APy'qyVly. Here 5sx is 
the Kronecker delta and e = 1 (^f: 1) gives the isotropic 
(anisotropic) Rashba-type interaction. We take the mo¬ 
mentum state to be a distribution in the x direction with 
e^Py p/s and the vacuum state in the y direction. With 
the internal state |p) = |a| |t) + \b\ ||), where |a| ^ |6| 

and (j) = "K12, such that {aP) = 0, we obtain 

(\ ePtP 

{x{t)) = [ ^ - 1 j T - [sin {2Pxt)] , (13a) 

{yij)) = ^ [cos {2Pxt) - 1]. (13b) 

The above equations lead to all the trochoids in Ref. 
[l3 - for appropriate parameters. When considering 
the Dresselhaus-type interaction ([6]) we obtain exactly 
the same time dependence as in Eq. m but with 
(ctx) Px+£ (cTy) py changed by {ax) py+£ {ay) px, affecting 
only the amplitude of the curves generated by Eq. (HU. 


In the second case we get a uniform motion in the s direc¬ 
tion and the expected trembling motion in the r direction: 

(s(r)) = (s(0)) -b (-1)1+^”" {as) T, (16a) 

(f(r)) = (f( 0 )) -b (- 1 )^+'^^“' r 

- Kr>c~^ {ps) {py + px) [sin {xt) - t] . (16b) 

Although Eq. (Il 6 bl) does not lead to trajectories similar 
to those discussed in Ref. [13 , we note that Eq. 0 can 
have its parameters properly adjusted so as to produce 
cycloidal motion.. 

Ionic Lissajous trajectories. Now we revisit the 
anisotropic version of the Rashba-type Hamiltonian (| 6 |). 
H = {'jxPxO'x — lyPy'^y)- By adjusting the laser fields 
such that » 7 ^ and preparing the electronic state as 
an eigenstate of cr^, we obtain 

(s(r)) = (s(0)) - {pr) [1 - cos {wt)] , (17a) 

(f(r)) = (f( 0 )) -b (psp/^) [1 - cos {Vwt)] , (17b) 
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FIG. 2: (a) Lissajou curves derived from Eq. (O with VO = 1.96 and riQ,zv and around 50 and 35 KHz, respectively, 

(b) Lissajous curves numerically computed from Eq. (IIOII . starting from initial coherent state 0 = 1 and electronic levels in the 
superposition \ip) = 14) — i IJ-); with N — 1 and 7 s/As = 1. In (c) and (d) we consider the same parameters as in (b), except 
for (c) 0 = 1 — i and (d) N = 2. In (e) we consider the same initial states as in (b) but N = 2, yaj/Aj, = 0.4 and 'yyjAy = 1. 
Finally, (f) follows from the same parameters as in (b) but considering dissipative mechanisms in the vibrational degrees of 
freedom, with a damping rate = 10 ~"^ 7 s/As. 


with w = {jx/jy) (Py)- In Fig- 2(a) we show Lissajous 
curves governed by Eq. (HZl) with zu = 1.96 and frequen¬ 
cies r]flw and rj^y/waioxaiA 50 and 35 KHz, respectively. 
In this case, the Rashba energy jy/px + Px is around 
10“^° eV. In what follows we obtain Lissajous figures 
from our previously introduced bounded SO interaction 
in Eq. (nni), which requires only the preparation of an 
initial vibrational coherent state |0) (differently from all 
the equations of motion derived above, which rely on the 
preparation of the initial vibrational state |'!/'(0))). In 
Fig. 2(b) the mean value of the ionic position has been 
numerically computed (running in QuTiP [l5l |l from Eq. 

with = 1, 7 s/As = 1, and starting from the 
vibrational mode in the coherent state 0 = 1 and elec¬ 
tronic levels in the superposition \ip) = |t) —* li)- In Figs. 
2(c) and 2(d) we consider the same parameters as in Fig. 
2(b) except for 0 = 1 — i in Fig. 2(c) and iV = 2 in 
Fig. 2(d). In Fig. 2(e) we consider iV = 2, "fxj^x = 0.4, 
"fy != 1, and the same initial states as in Fig. 2(b). 

Detrimental effects. To show that our calculated 
ionic trajectories are robust, we have included damp¬ 
ing effects due to the environment. Figure 2(f) 
is similar (same parameters) to in Fig.2(b) but ac¬ 
counts for dissipative mechanisms in the vibrational 
degrees of freedom as described by the Lindbladian 


(C.s/2) [2a.,/oot — ata.sP — pala.y ], with a damping rate 
Cs = 10“^7 s/As [l^. Clearly, the trajectories are robust 
and visible for realistic parameters. 

We have presented a protocol for generating Lissajous 
curves with the vibrational motion of an ion in a two- 
dimensional trap. It relies on the unique capability of our 
setup to realize Rashba- and Dresselhaus-type SO inter¬ 
actions, which allows us to simulate solid-state SO effects 
within a highly controllable trapped-ion experiment. We 
have also verified that Lissajous curves can be derived 
from upper-bound SO interactions, which may bring new 
perspectives to the subject. Addressing some interesting 
issues to be investigated further, we first observe that a 
straightforward extension to the case of many trapped 
ions, where the strong and tunable (up to W^Hz) SO 
stren gth can be used to explore quantum phase transi¬ 
tions [13 and chaotic behavior [18|. Finally, motivated 
by the results above, we believe it is worth to investigate 
the role of the “Bounded” spin orbit interaction in solid 
states systems [l9| . 
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